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1 Introduction 

In recent years the exploration of the universe has provided detailed information about 
the distance-redshift relationship of many sources at very small redshifts as well as at 
redshifts of order unity. These observations, if interpreted within the framework of 
a homogeneous Friedmann-Lemaitre-Robertson- Walker (FLRW) metric indicate that 
the universe is presently undergoing a phase of accelerated expansion [1]. Such an 
accelerated expansion is most commonly interpreted as evidence for the presence of 
a negative pressure component ("Dark Energy") in the mass-energy density of the 
universe. This has given rise to the so-called "concordance" fiat AC DM model, where 
dark energy constitutes 70% of the present energy budget in the form of a cosmological 
constant along with 30% of matter. A second fact lending credence to the AC DM 
model is the observation that CMB excludes, with high precision, the presence of 
spatial curvature. Applying this to the late time universe in an FLRW model leads 
to a mismatch between the observed matter density and the FLRW critical density. 
This mismatch can also be cured by the presence of a Dark Energy component. 

However, from the point of view of fundamental physics the presence of such a tiny 
cosmological constant is almost absurd: first, the scale is extremely low compared to 
any possible fundamental scale, and second it is tuned in such a way that it shows 
up just at 2; ~ 1, where we live. This is so baffling'' that before concluding that we 
are indeed facing this huge puzzle, it is worth making efforts to explore if the correct 
equations are being used to fit the data. 

The homogeneous FLRW metric certainly offers the simplest paradigm to interpret 
cosmological observations. Indeed, the FLRW model is a very good approximation 
in the Early Universe as probed by the homogeneity of the CMB (density contrasts 
in photons and dark matter are of the order 10~^ and 10~^ respectively). However at 
low redshifts, the density contrast in matter grows up to values of 0{1) and beyond: 

"^There have been some ideas/efforts to address these issues. For instance, see string-gas inspired 
coupled quintessence models [2] by one of the authors. However, it is fair to say that there is not 
yet a completely satisfactory way of solving both the "smallness" and the "coincidence" problems. 
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small scales become nonlinear first, and then more and more scales enter this regime. 
Today, scales of order of 60//iMpc (which is 2% of the size of the horizon) have an 
average density contrast of order 1. Structures observed with density contrast larger 
than 0.1 (so, already in a mildly nonlinear regime) extend to few hundreds of Mpc 
(around 10% of the size of the horizon). It is thus not clear or obvious that one can 
keep using FLRW metric to interpret the large amount of high precision experimental 
data collected in recent times in a regime where most of the mass in the Universe is 
clumped into structures or it is forming structures. 

Three are three main physical effects that we could be missing while using the 
FLRW model: 

I. The "overall" dynamics of a Universe with these inhomogeneities could be sig- 

nificantly different from the FLRW dynamics. 

II. Even if the dynamics is approximatively the same, one has to wonder about the 

light propagation in a clumpy Universe: since all our conclusions about the 
dynamics are based on observations of light, this is of crucial importance. 

III. The fact that we have just one observer could influence significantly the obser- 

vations: we could live in an especially underdense or overdense region, which 
may induce large corrections on the observations. 

In recent times there has been an on-going debate as to whether some of these effects 
may be large enough to give rise to an "apparent acceleration" leading us to believe 
in dark energy when there is none. The aim of this paper is to construct an exact 
model of non-linear structure formation where we can systematically study these 
effects. We wanted a model where one can study (possibly also analytically) the 
effects that inhomogeneities can cause on quantities like the redshift and the angular 
distance, thereby hopefully providing us with valuable insights to resolve the issue. 
The full problem of solving Einstein equations and light propagation for a realistic 
mass distribution is beyond present human possibilities, so some simplification is 
needed. Our simplification is to assume spherical symmetry and in particular our 
model is based on Lemaitre-Tolman-Bondi (LTB) metrics which are exact spherically 
symmetric solutions of general relativity with only dust. This approach has the 
obvious advantage of being exact: it does not use perturbation theory (as most of 
previous literature did to address I), it allows to study light propagation without 
any FLRW assumptions (as opposed to literature addressing point II), studying the 
corrections on the distance and also on the redshift of light (which has been ignored 
by most previous literature). The challenge however, is to make these kind of models 
as realistic as possible, or at least be able to extrapolate from them the information 
that is relevant for our Universe. 

So far such models have been used mostly to understand what happens to obser- 
vations if we are living in a special position, a void for example [3, 4, 5, 6, 7, 8, 9, 
10, 11, 12, 13], and therefore can at best be thought of as describing the local uni- 
verse. It is interesting to go beyond this approach. Firstly, we may be missing some 
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important effects (such as light propagating through several structures) by looking 
just at the local universe. Secondly, although these models show interesting effects 
capable of mimicking dark energy, it is not clear how realistic are these setups [15, 16] 
(also, some of these results have been questioned by [17]). On the other hand, we 
note that a different approach (aimed at solving question I), based on back reaction 
of inhomogeneities on the overall dynamics of the Universe, suggests a much smaller 
correction at the first nontrivial order [18, 19, 20, 21], although even in this frame- 
work some authors have argued in favour of a much larger effect [19, 22, 23]. Thus 
it is imperative that we understand the origin of the large effects and in the process 
hopefully also bridge the gap between the perturbative back-reaction effects and the 
exact toy model approach. 

It is clear that to address these issues we have to go beyond a "local" or a "per- 
turbative" description. Thus, although we work with the simplifying assumption of 
spherical symmetry, we try to make our model realistic in the following sense. The 
universe looks like an onion (and we frequently refer to it as the "Onion model"): 
there is a homogeneous background density and on top of it density fluctuations as a 
function of the radial coordinate. Our model preserves large-scale homogeneity and 
the shells are also distributed in a homogeneous way, making the picture a lot closer 
to the real world. One may wonder, nonetheless, that the centre is still a special 
point. In fact, we find that the centre has some special properties. For instance 
typically, the curvature tends to become large at the centre. In this case if we find 
significant corrections to, say, the redshift-luminosity distance relation near the cen- 
tre, it is hard to distinguish whether the correction originates from the non-linear 
structure formation, or the "excess curvature" . This therefore can jeopardise what 
we are trying to avoid, that of violating the cosmological principle. We deal with 
this problem by considering an observer who sits in a generic position and looks at 
sources in the radial direction, along which we have periodic inhomogeneities. To our 
knowledge, the luminosity-redshift relation in LTB models has always been analyzed 
putting the observer at the centre of the coordinates^. For the first time we derive 
an exact expression for the luminosity distance of an object as seen by an off-centre 
observer^ . 

At early times (at last scattering for example) we assume the density fluctuations 
to have an amplitude of the order of the fluctuations observed in the CMB. Then 
we are able to follow exactly the evolution of density fluctuations. At redshifts of 
z > 0(10) the density contrast grows exactly as in a perturbed FLRW. Then, when 
the density contrast becomes of 0(\) nonlinear clustering appears. So at low redshift 
the density does what is expected: overdense regions start contracting and they 
become thin shells (mimicking structures), while underdense regions become larger 
(mimicking voids) and eventually they occupy most of the volume. One can see 
this completely analytically upto high density contrasts (we have checked that the 
analytical solution is good even with bpjp ~ (9(100)), which already makes it an 

^Except for an analysis at small z in [24]. 

^In this respect we note that [17] has shown that locally around the centre in the LTB model the 
deceleration parameter cannot be negative. However, this result is only local and so the universe 
may accelerate slightly far away from the observer. 
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interesting model for the theory of structure formation. However, the really nice 
feature of the Onion model is that one can even calculate how the redshift and the 
various distances get corrected when we are in this non-linear regime. 

In this model we can solve for the radial geodesies, therefore determining exactly 
what happens to the redshift of light. We place the observer in a generic location and 
measure the redshift as a function of the source coordinate. Then we trace a beam 
of light that starts from the source with some initial angle, and in particular we are 
able to compute the area seen by the observer. This provides us with the angular 
distance (-Da) or equivalently the luminosity distance (-Dl) for a generic source. In 
this way we provide the — z relationship (or the Da — z relationship) exactly. 
Thus we are able to study the three possible effects (I, II, and III) both numerically 
and analytically. We now discuss our findings briefly. 

The first issue (I) of whether there is an overall effect in expansion due to inho- 
mogeneities has been investigated only recently. Buchert [25] has defined an average 
volume expansion and has developed a formalism to take into account of the back- 
reaction of the inhomogeneities, although the problem of computing the size and 
the effect of the corrections is left unanswered. Rasanen [19] has tried to estimate 
the effect of inhomogeneities (extending the formalism developed in [26], albeit in a 
different context) using perturbation theory, finding an effect of order 10"^ (in agree- 
ment with [18]) and speculating about nonlinear effects. A complete second order 
calculation has been done in [20], and in a subsequent work one of us has shown [22] 
that the perturbative series is likely to diverge at small z due to the fact that all the 
perturbative corrections become of the same size (each one is 10~^), and one may 
hope that these corrections give rise to acceleration (which it is shown to be in prin- 
ciple possible [27, 22]). The fact that this may lead to acceleration has been explored 
in subsequent papers [23] and criticized in others ([21] for example have argued the 
effect to be too small). Unfortunately a perturbative approach seem inadequate to 
answer this question as the phenomenon is intrinsically nonlinear. 

On the other hand, since our model is an exact solution of Einstein equations, we 
can expect to get faithful results. Based on our analysis we do not find any significant 
overall effect^. Quantities, such as the matter density, on an average, still behave as 
in the homogeneous Einstein-de Sitter (EdS) universe {i.e. a matter dominated fiat 
Universe) . We remind the reader that if the universe starts to accelerate the matter 
density should dilute faster than EdS models. 

The second issue (II) has been investigated in the past already in the '60s. The 
main point here is that the light that travels to us in a realistic Universe is not red- 
shifted in the same way as in a FLRW metric, nor the angles and luminosities evolve 
in the same manner: it is more likely that a photon travels through almost empty 
space rather than meeting structures which are clumped in small sizes. However, 
typically work in the past has focused on what happens to the angles (or equivalently 
on luminosities), almost ignoring the problem of the redshift. Early pioneering work 
on this subject was due to Zeldovich [28] that estimated the effect on the angular 

^In agreement with the similar findings of [15, 48] 
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distance (Da) on light travelling through an empty cone embedded in an FLRW 
metric. A similar treatment based on [29] has been used recently to analyze data 
by the Supernova Cosmology Project Collaboration (see [30], fig.8). While giving a 
quite relevant correction (for large matter fraction Qm), the effect does not allow to 
get rid of the cosmo logical constant. Therefore the collaboration [30] estimates to 
have control on the inhomogeneities. However such calculations are incomplete: the 
dynamics is completely ignored, and more importantly the redshift experienced by a 
light ray is assumed to be the usual FLRW result. 

Instead in the Onion model we could keep track of the corrections to both the 
redshift and the luminosity distance and in fact we found both of them to be signif- 
icant, once the inhomogeneities become large. We found that the correction in the 
redshift with respect to a homogeneous model is of the same order, and even larger, 
than the correction in luminosity distance. Qualitatively one can understand where 
such corrections come from: in the underdense regions the expansion is faster and 
the photons suffer a larger redshift than it would in the corresponding EdS model, 
while in the overdensities the expansion is balanced by the gravitational collapse and 
the photons experience milder redshift (or even a blue shift!). Here comes however 
one of the shortcomings of the Onion model: a radial light ray in the Onion model 
unavoidably meets underdense and overdense structures and they tend to average out 
the corrections. In the real Universe, instead, where the light hardly encounters any 
structure, one does not expect such cancellations to occur. Nonetheless, based on 
these considerations, we can try to go beyond the Onion model: since the photon in 
the real world is mostly passing through voids it gets redshifted faster as the non- 
linearities increase with time, thereby potentially producing the effect of apparent 
acceleration. In this paper we only present an estimate of such an effect leaving a 
more detailed study for future work. 

Next we come to the third issue (III) of whether the position of the observer 
can have any important consequences. In general, in our analysis, we found that the 
corrections to redshift or luminosity distance is controlled by two quantities (in fact 
it is the product of the two): the amplitude of density fluctuations, Sp/p and the 
ratio L/rhor, where L is the wavelength of density fluctuations and rjior is the horizon 
distance. As we mentioned earlier, the length scale at which the density fluctuations 
become non-linear, i.e. Sp/p ~ the latter ratio is only around a few percent. 

Therefore naively one may expect these corrections to be too small to be relevant 
for supernova cosmology (as argued in several literature). However the issue is more 
subtle: near the observer when the redshift and distances are themselves small these 
corrections become significant as they do not proportionately decrease. Instead of 
the ratio L/r^or, the relative correction is now governed by the ratio L/6r, where 
Sr is just the distance between the source and the observer. Thus in the first few 
oscillations, especially the first one, when (5r ~ L we obtain significant corrections. 

In our analysis we find that, if we want to explain the Hubble diagram solely by the 
effect coming from the first oscillation, then the location of the observer is of crucial 
importance^. Basically, what one needs to explain is the mismatch in the measurement 

*This situation can change drastically if we are also able to account for the second effect, that of 
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of the local Hubble parameter (between redshifts 0.03 < z < 0.07) and the luminosity 
distance observed for high redshift supernovae (typically between redshifts 0.4 < z < 
1.5). Although we can place the observer in any generic location, we find that the 
best choice, to account for the larger locally measured Hubble parameter, is if we 
place the observer in an underdense region. Our main result is that we found that 
we require 3 — 4 times larger local density contrast than the average value inferred by 
CMB measurements. With these assumptions, however, we find that the Onion model 
can reasonably fit the supernova data and in this sense our model corroborates some 
of the findings of the models based on local voids. We further find that the Onion 
model can be consistent with other observations such as local density measurements, 
the first peak position of the CMB (which measures the curvature of the universe) 
and baryon acoustic oscillations. 

The Onion model also offers us the opportunity to look at other relevant observa- 
tions. For instance, since we have the freedom to choose the position of the observer 
we can in principle calculate quantities such as anisotropies in the measurement of 
Hubble constant which are essential to make a realistic assessment of the viability of 
any model relying on "us" being in a special position. 

Finally, suppose that inhomogeneities do not give a large enough effect to mimick 
Dark Energy, we can still estimate what is a minimal effect due to structures that 
typically affects distance measurements in the real Universe, i.e. an universe con- 
structed by just propagating the CMB fluctuations to today. We find that a typical 
effect gives a correction of order 0.15 apparent magnitudes to the Hubble diagram 
at all redshifts. We can also explain another source of uncertainty in supernovae: 
the observed nearby supernovae seem to have larger intrinsic scatter [31] which can 
easily be explained by the large oscillations that we get in the Dl{z) curve close to 
the observer. We briefly discuss such applications in this paper but mostly focus on 
the problem of dark energy. 

The structure of the paper is as follows. In section 2 we introduce our LTB model, 
and in section 3 we solve it and discuss structure formation. In section 4 we solve for 
light propagation in this background, both analytically and numerically. In section 5 
we discuss how to find the luminosity distance, both analytically and numerically. In 
section 6 we first construct a setup that could mimick Dark Energy, and then ask how 
realistic is such a setup, discussing observations such as local measurements of matter 
densities, CMB acoustic peaks and baryon oscillations. We then go on to discuss light 
propagation in the real Universe as opposed to the Onion model, and we estimate 
how large can the missing effect be. We also estimate minimal uncertainties that 
inhomogeneities can lead to in the Dl — z relationship. Finally in section 7 we draw 
our conclusions. Appendices contain many technical details: Appendix A contains a 
discussion of the analytical approximations that we employ. Appendix B contains a 
brief discussion of what happens near the centre of the LTB model. Appendix C and D 
contain the analytical derivation of (respectively) the time and the redshift along the 
photon trajectory and Appendix E contains the full calculation of -Dl for an off-centre 
observer. In appendix F we show how our results for the off-centre angular distance is 

light mostly propagating through voids. 
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consistent with the optical scalar evolution equations for the angular distance derived 
by Sachs in 1961 [32]. 



2 The model: LTB Universe 

As we mentioned before, our model is based on an LTB metric. LTB metrics are 
spherically symmetric solutions of Einstein's equations with only dust (pressureless 
matter) acting as the source of energy density. So, in this respect, they are similar to 
the EdS model, but one is allowed to have an inhomogeneous distribution of matter 
along the radial direction, and this is what we are going to exploit in order to describe 
structure formation. Such a solution was first proposed by Lemaitre and it was later 
discussed by Tolman and Bondi [33]. Early applications of LTB [33] models were 
concentrated on the study of voids [34, 35]. Raine and Thomas used a LTB model 
to study the CMB anisotropics produced by an overdense matter distribution having 
a small density contrast of about 1.5 percent on a scale of order 1000 Mpc [36]. 
This work has then been generalized to the non-linear cases in [37]. The problem of 
light propagation in an LTB model and its impact on the Hubble and deceleration 
parameter was dealt with in [4, 38]. 

After the release of the new data on SNIa, following Celerier's [3] interesting 
paper exploring the possibility that large scale inhomogeneities may be able to mim- 
ick dark energy, several authors have used LTB solutions to understand such ef- 
fects [6, 7, 8, 9, 10, 11, 12, 13, 14, 17]. Some authors have tried to study the effect 
of averaging the Einstein equations for LTB inhomogeneous solutions and its result- 
ing backreaction which may also mimick an accelerating universe [39, 12, 13, 40], 
while [15] has questioned such approaches, as requiring unrealistic density fluctua- 
tions. Our Onion model however differs from these "void models" in the following 
sense. Firstly, our LTB solution incorporates the entire Universe, without giving any 
physical meaning to the existence of a centre. Secondly, the philosophy behind our 
model is very different as compared to the Void models. We start with an initial 
universe which resembles the CMB epoch and then evolve to our present era, while 
the void models are typically valid only at late times. Also, very importantly we are 
able to analytically understand all our results thereby clarifying several issues that 
have been debated in the literature. 

There are, however, limitations to our model of structure formation. The most 
apparent limitation is that we can describe only 2-dimensional structures (that look 
like thin spherical shells). While this is not completely realistic, it is nonetheless 
similar to what happens in the study of structure formation: using for example the 
so-called Zeldovich approximation [41] it can be shown that structures form after the 
collapse along one of the dimensions, leading to two-dimensional pancake-like objects. 

The really crucial limitation of our model instead is that, as in the Zeldovich 
approximation once these structures start collapsing, there is nothing that can prevent 
them to become infinitely thin (shell crossing). This is a serious limitation since, in the 
real world, objects which reach high density at some point can virialize and form stable 
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structures. The reason why our model does not include this feature is that we have 
motion only in the radial direction (and no pressure whatsoever) and therefore there 
is nothing that can prevent collapse: in the real world rotational motion is allowed, 
and systems can virialize if centrifugal forces balance the gravitational attraction. 

In order to avoid this singular collapse, we choose a scale with an initial density 
contrast such that the collapse happens only in the future with respect to observation 
time. This forces the redshift, Zni, at which the density become nonlinear, to be quite 
close to zero. In other words we cannot track non-linearities for too long time (upto 
around 2; ~ 1) 

We will come back to discuss these limitations later in section (6.2), but for now 
our aim is to construct an LTB model which can resemble a realistic universe with a 
largely homogeneous component and evolving density fluctuations on top of it. 



2.1 LTB metrics 

We start with a brief review of LTB solutions. The LTB metric (in units c = 1) can 
be written as 

ds^ = -dt^ + S\r,t)dr^ + R\r,t){d9^ + sm^9dip''), (2.1) 

where we employ comoving coordinates (r, 6, ip) and proper time t. 

Einstein's equations with the stress-energy tensor of a dust, imply the following 
constraints: 

1 -P, , GM(r) , , , , 

M'(r) 

where a dot denotes partial differentiation with respect to t and a prime with respect 
to r. p{r,t) is the energy density of the matter, and G = l/m?pi is the Newton 
constant. The functions E{r) and M(r) are left arbitrary. The function E{r) is 
related to the spatial curvature, while M(r) approximately corresponds to the mass 
inside a sphere of comoving radial coordinate r [3]. 

One easily verifies that (2.3) has the following solutions for R{r,t), which differ 
owing to the sign of the function E{r): 
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• For E{r) > 0, 



R 

t - h{r) 



■1), 



(2.5) 



E{r) = 0, 



R{r,t) 



'9GM{r 



n 1/3 



[t-t,(r)]3 



and E(r) < 0, 



R 

t - U{r) 



GM{r) 
-2E{ry 
GM{r 

i-mr)] 



1 — cosm) 



(2.6) 



3/2 



— sm u) , 



where tfe(r), often known as the "bang-time", is another arbitrary function of r. It 
is interpreted, for cosmological purposes, as a Big-Bang singularity surface at which 
i?(r, t) = 0. This is analogous to the scale-factor vanishing at the big bang singularity 
in the homogeneous FLRW models. One can choose tfe(r) = at the symmetry centre 
(r = 0) by an appropriate translation of the t = constant surfaces and describe the 
universe by the t > tb(r) part of the (r, t) plane, increasing t corresponding to the 
future direction. 

To summarize, LTB models contain three arbitrary functions, M(r), E{r) and 
th{r). It is easy and instructive to check how to recover the homogeneous limit. First 
of all, one has to choose the big bang time to be the same everywhere: conventionally 
tb{r) = 0. Then the E{r) function has to be chosen as: 

E{r) =-Ey Closed FLRW 

eIt) = Flat FLRW (2.7) 

E{r) = EqP Open FLRW 

where Eq is any positive number. One immediately recovers the homogeneous FLRW 
models, with Eq parameterizing the amount of curvature in the model. 



2.2 Choice of M(r), E{r) and tb{r) 

First of all, the reader should not be confused by the terminology which is usual in 
the literature: models with E{r) < 0, E{r) = 0, E{r) > are called respectively 
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open, flat and closed LTB. However as we are going to show, a model with E{r) > 
can mimick easily a flat FLRW model, and can be consistent with measurements of 
the flrst CMB peak that tell us that "the Universe is flat" . 

As a flrst remark we note that one of the three functions that describe an LTB 
model is just a "gauge" degree of freedom. In fact, by an appropriate rescaling of the 
radial coordinate r the function M(r) can always be set to any monotonic function. 
For example, starting with a radial coordinate r and a generic function M(f), one 
can perform a coordinate transformation r = r(r) without affecting the form of the 
metric, i.e. 

M{r) = M(?), E{r) = E(j) and i?(r, t) = R(r, t) . (2.8) 
Thus we are free to deflne a new coordinate as follows: 

AMl)y"^MM = i^Af„V. (2.9) 



r r 



where Mq is an arbitrary mass scale. In the rest of the paper we always use this 
convention. Thus, we are left just with two functions: E{r) and tb{r). 

It is easy to see that the function th{r) becomes negligible at late times: t S> 
tb{r). Therefore the function tb{r) can only represent a decaying mode in the density 
fluctuation^. Since the information contained in tb{r) gets lost at late times, we can 
just ignore this function and set tb{r) = for all r. Instead, in order to describe mass 
condensation and void formation at late times, we have to choose a non-trivial E{r). 

At this point we have the choice between E{r) > or E{r) < 0^^. It is also 
possible in principle to have a mixed model with some regions where E{r) < and 
other regions where E{r) > 0. For simplicity, we make the choice to always have 
E{r) > 0. As we will show, this can describe mass condensation and at the same 
time be consistent with the CMB flrst peak measurements. 

For later convenience we enumerate the simplifled equations for our model with 

tb{r) = 

R{r, t) = r ^ > (cosh^- 1) , (2.10) 



where we have deflned 



M = ^. (2.12) 
mpi 



^It is easy to check that no matter what function we use as ti,{r) the universe becomes more and 
more homogeneous with time (with trivial E{r)), which is exactly the opposite of what we want. 
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Choosing E{r) ~ gives us back the flat homogeneous EdS model. 
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3 Structure Formation 



In this section we will try to understand whether open LTB metrics can mimick ho- 
mogeneity over large distance scales, as well as describe inhomogeneities (representing 
galaxies, clusters of galaxies, etc..) appropriately. In other words we will try to single 
out the i?(r)'s which can "approximate" the density profile of our universe. 



3.1 Large scale homogeneity and "small m" approximation 

Throughout the entire paper we will be employing what we call the "small u" ap- 
proximation. In appendix A we show how this is a very good approximation for a 
space-time region large enough to encompass our entire causal horizon. The strength 
of the approximation lies in the fact that it can accurately describe the dynamics 
even when Sp/p 3> 1, and thereby it lets us study the effect of non-linear structure 
formation on the different physical quantities. In fact, it is mainly the validity of this 
approximation along with an adiabatic approximation scheme (see appendix C for 
details) that provides us with so much analytic control. 

In (2.11), if we only keep next to leading terms in u we have 



and 

6E{r) V 2 4! 

It is clear from (3.1) that this approximation is valid {i.e. u is small) when 

V 27rM2r3 J Mr \ tt J 

Now, if we choose an E{r) that grows less than quadratically then, as long as we 
are sufficiently far away from the centre, we can trust the small u approximation 
(see appendix A for a more precise estimate). We will see later that, in fact, we do 
not want an E{r) growing faster (or equal) than because it spoils the large-scale 
homogeneity or it introduces unwanted overall spatial curvature. 

Proceeding accordingly, (3.1) yields 

,2 



which can be solved giving us 



60 



y3 



u-v--. (3.4) 
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Combining (3.2), (3.3) and (3.4) we find 

R{r,t) ^ (67r)^/V(Mt)i 
where we have defined 



1 + a 



E{r){Mt)^ 
{rMf 



a = n 

V 40007r2 



1/3 



(3.5) 



(3.6) 



As one can easily check, the prefactor in the expression for i?(r, t) corresponds to the 
flat FLRW result while E{r) corresponds to corrections in R{r,t). 

Next, let us try to understand how the density profile looks like. First we will 
approximately compute the average density inside any sphere of radius f. 



(P) 



M,o, _ / drR^R'p{r,t) 



(3.7) 



Kot / drR'^R' ° / dRR^ 

where we have assumed that E{r) <^ 1 (see appendix A for a justification). The last 
equation can be written as: 

(p) = — i + a 



67i{Mtf 



E{f){Mt) 



{fMf 



(3.8) 



Imposing large scale homogeneity basically means that (p) should have a well 
defined limit as f — oo. From the above expression it is clear that this happens if 

E(f) 



lim 

r— >oo 



[fMy 



const. 



which means that E{r) cannot increase faster than as r — ^> oo. Moreover, as we 
have pointed out in the previous section, a quadratic E{r) corresponds to a model 
with constant spatial curvature (2.7). This is strongly constrained by CMB data [43], 
because the angular diameter distance (-Da) for an universe with spatial curvature at 
high redshift is very different from the flat case prediction, and the measurement of 
the first peak position strongly favours the flat model (or more precisely constrains 
the amount of constant spatial curvature to be very small). Therefore, in order to 
avoid curvature, E{r) actually needs to increase slower than as r — > oo. We will 
now further show that only a linearly growing E{r) (times a factor that includes fluc- 
tuations) guarantees a realistic density profile (relative heights of peaks and troughs of 
the fluctuations) that does not change as we move outwards from the centre towards 
infinity. We will explicitly see in section 6.1, that such an E{r) is also compatible 
with having a flat universe. 



3.2 Inhomogeneities 

To compute p{r,t) we have to compute R' . From (3.5) we find 

'E(r 



R'{r,t) = {QTif'\Mty^ 



1 + a(Mt)3 



rM2 



(3.9) 
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From (2.4) we then find 



p{r, t) 



1 + a 



E(r){Mt)^ 



'[l+a(M«)ifai)' 



(3.10) 



We note that the correction term in the second factor of the denominator is actually 
proportional to v"^ 



a- 



E{r){Mty. 
{rMf 



a 



(9V2) 



2/3 



v' < 1 



and therefore can be ignored consistently with our "small m" approximation. Thus 
we have 



p(r, t) 



67r(Mt)2 l + a{Mt)3A{r) 



where A{r) 



E{r) 



(3.11) 



A few comments are in order. First, observe that the FLRW behaviour for the 
density is given by the prefactor multiplying the denominator. The fluctuations are 
provided by the presence of A{r) in the denominator. It is also clear that if we do 
not want the relative peak to trough ratio to change, then the function A{r) must be 
bounded 

A-min ^ A(r) < Amax ■ (3-12) 

In other words there should not be any "overall" growth of A{r) as r increases. From 
the definition of A{r) (3.11) this means that E{r) has to have an "overall" linear 
growth apart from the oscillatory fluctuations; if E{r) grows faster or slower than 
linear, A{r) starts to grow or shrink respectively. 

It is clear that the density contrast (5) defined in the usual manner 
is controlled by 

e{r,t) = a{Mt)lA{r) . (3.14) 
One obtains a simple relation, ignoring the small correction of (3.8) 

S = _^ = _^^mi^. (3.15) 

1 + e l + a{Mt)lA{r) 

Let us assume, for the purpose of discussing the dynamics, that Amm is negative. We 
will see in the next subsection that this will be required in order to avoid curvature. 
It is clear that the density fluctuates between 

a{Mt)lA^a, ^ ^ ^ a(Mt)i(-A^,„) _ ^ 

Omin{t) = ~ . < < ^ 2 = Omax{t) ■ (3.16) 

1 + a{Mt)sA^ax 1 + a(Mt)3 A™„ 
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When A{r) is close to Amax we have a void, while when it is close to Amin, it signals 
an overdensity. It is nice to see how the transition from linear to nonlinear regime 
occurs. When e -C 1, 5 ~ e and we are in the linear regime: an initial density fluc- 
tuation grows at small times as ts, in agreement with the prediction of cosmological 
perturbation theory. When e is no longer small the density contrast grows rapidly, 
and it is comforting to see that this result is the same as found within the Zeldovich 
approximation [41]. 

Even if one chooses = —A^ax, in the nonlinear regime typically |5mm(^)| < 
Smax{t)- This is because the denominator in (3.15) comes in play. For positive e 
the denominator does not change much, but for negative e as it approaches —1, the 
denominator tends to zero thereby giving rise to a very large 6max- In fact, as t 
increases there comes a point when the denominator in (3.15) precisely vanish so 
that the density in the over-dense region is infinite (Smaxit) oo), while Sminit) 
remains finite. The "critical" time when this happens is given by 



In the real Universe, rather than the density going to infinity, the system typically 
virializes, forms rotating structures (galaxies, clusters, etc.), and falls out of the 
Hubble expansion. As a result, Smax{t) in reality should saturate to a maximal value 
(unless a black hole is formed), while the underdense regions keep forming bigger 
and emptier voids. Unfortunately in our model we cannot incorporate the effects 
of virialization and cannot trust the results beyond tg. We will come back to the 
important question of whether this limitation may lead us to underestimate the effects 
of inhomogeneities on the luminosity-redshift relation in section 6.2. 



Let us come back to the main issue of finding appropriate £'(r)'s that can describe 
the large scale structures. We have seen how the function A{r) controls the density 
profile and in particular that we need a bounded A{r). From now on we will focus on 
a model with sinusoidal density contrasts. (It is possible to consider a spectrum of 
many wavelengths, but in order to keep the analysis simple we will focus on a single 
wavelength L.) As we will justify later, to achieve this it is sufficient to obtain a 
density profile that is sinusoidal in the comoving coordinate r. 

It is evident from (3.14) that, to obtain such periodic oscillations of over and under 
dense regions, we require a periodic A{r): 



(3.17) 



ts = 




3.3 The "Onion" Model 





Integrating (3.18) we find 



E{r) = r'^Ao + tAq - r 



AiL 



cos 




(3.19) 
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The term proportional to Aq gives rise to curvature and hence we do not consider it 
anymore. It is easy to see that this is precisely the term that one would generate if 
Amin > and this is why we said before that we needed Amin < 0. In the absence of 
^0, E{r) becomes 



E(r) 



rAa — r cos 

271 



27rr 



We note that the term proportional to Aq does not enter into the density profile and 
therefore we do not expect it to play any significant role in modifying the luminosity 
distance-redshift relation (we discuss this in sect. 4.1). Since we want E{r) to be 
always positive, it seems convenient to choose 



so that 



E(r) = M'r 



2 



2ti 



An 



1 — COS 



All 
' ~2^ 

2'Kr 



L 



~ 2 . ^2 fur 

M^r sm^ — 

271 \ L 



(3.20) 



(3.21) 



This defines our model uniquely, with two parameters, L and Ai. For later conve- 
nience we enumerate R and R' for our Onion model: 



R{r,t) ^ (67r)^/V(Mt)t 



e L 9 /vrr\ 

1 + sin^ ( — ) 

271 r \L J 



and 



1 + e sin 



R'ir,t) = {Q7if/\Mt)-^ 
where we have introduced the "amplitude of oscillations 

e(t) = aAi{Mt)^ . 



271 r 



L 



(3.22) 



(3.23) 



(3.24) 



It is useful to observe that at any given time-slice, R is always proportional to 
the coordinate distance r, especially at earlier times when the fluctuations due to 
inhomogeneities are small. Now, it is clear from the metric that R is actually related 
to the proper distance; it is in fact exactly the space-like distance between two points 
along the same radial ray in a time slice. Thus to a good approximation, barring a 
proportionality factor, one can indeed treat the coordinate distance r as the "physical 
distance". This explains why it physically makes sense to consider a density profile 
which is just periodic in r. 

For an E[r) given by (3.21) the density profile looks like 

P ^, t) = F 7^r^ . 3.25 

'^^ ' 67r(Mt)2 [i + e(t)sin(^)] ^ ' 

with periodic over-dense and void regions like an onion. The amplitude of density 
contrast is controlled by the function e(t) and increases as time progresses. 
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Figure 1: The density contrast for a given L but at different epochs as a function of 
the radial coordinate, numerically (black dashed line) and analytically (orange solid 
line). The two curves are practically indistinguishable. 



As seen in figure (1) our analytical estimate agrees with the exact numerical 
solution very well. As we discuss in details in appendix A, this is possible because we 
can find a sufficiently large region in r — t plane where both E{r) and u are small. 

As anticipated, we encounter a problem of singularities, i.e. there is a time, tg, at 
which the density contrast blows up. In our model 

= -Arnax = -Ai with > , (3.26) 

so the time ts is given by (3.17): 

'.M=(J^- (3-27) 

On the other hand we know the density contrast at CMB, as measured first by 
COBE[42, 43], is given by (3.15) 



2 



ScMB = eitcMB) = aAi{MtcMB)~^ ■ (3.28) 
We will abbreviate "CMB" by "C" from now on. Thus we find 

r = ^c'" ^ (3-29) 

which means that the "singularity" time does not depend on the scale of the inho- 
mogeneities (L), but only depends on the initial amplitude of density fluctuations. 
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Clearly, this is also the time when the density contrast is just relevant for structure 
formation 6p/ p ~ 0{1), so we will choose t<j to lie close to to-, the observation time. 

We stress again here the limitations of our model. Due to the fact that we work 
in the region t < ts, we can have structure formation only at very low redshifts. In 
this sense the only scales that we can describe realistically are those that collapse 
roughly at the present epoch, that is with L of the order of several Mpc/h (where 
h is related as usual to the measured value of the Hubble constant Hq by Hq = 
100 Kmsec~^Mpc~^), that corresponds to the scale of the largest seen structures 
like voids or superclusters. Moreover if we work at times t < tg, where e{t) < 1, it 
is clear from (3.15), that the minimal density contrast that can be achieved in the 
Onion model is only 5 = —0.5. 



4 Radial Photon Trajectories and Redshift 

In the previous section we have seen how our LTB Onion model can be used to describe 
the nonlinear large scale structure of our universe. Our aim now is to find corrections 
to the Dl — z curve coming from inhomogeneities like voids and superclusters that 
we observe in our present universe. 

In the Onion model structures are concentric shells of high density, and voids are 
replaced by the space between the shells. The physical setup that we consider is to 
have the observer far from the centre of the LTB patch in order to avoid as much 
as possible any peculiar behaviour due to the fact that the centre is a special point. 
Moreover in this way we can make the photon propagate in a region in which u is 
always small, and under good analytical control. So, the observer will be located 
at some nonzero tq, and we will choose it to be large enough such that u is small, 
but not too large so that E{r) also remains small. Then we will consider sources 
at variable position r > tq, but always along a radial trajectory. In this case we 
are able to easily solve the geodesic equations numerically and also to find analytic 
approximations. In the same configuration we will be able to define and compute the 
angular (and luminosity) distance in section 5. 

At this point one may wonder about other directions that an observer can look 
at. We do not consider them because approximately the picture looks the same and 
we do not expect any qualitatively new physics. It is another matter that technically 
it is extremely hard to solve for such non-radial trajectories. 

4.1 Null Geodesic Equations 

The first step towards obtaining a relationship between the luminosity distance and 
the redshift is to solve the null geodesic equations. In other words, we have to obtain 
t{r) for a photon trajectory converging at say r = ro at time t = to- (The subscript 
O from now on will denote the observer). 
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The equation for this radial photon trajectory is simply given by 



dt{r) R'{r,t{r)) 



dr y/1 + 2E(r) 



(4.1) 



First, note that we recover the FLRW limit by letting E{r) here and in (3.5). 
In this case we find 

t(r)V3 ^ ^1/3 _ ^^2/3^^ _ ^ ^l/3(^) ^ (4 2) 

where 

(670^3 
P - 3 • 

The nice feature of the Onion model is that it allows us to find an analytic approxima- 
tion to (4.1), in the small u and small E framework (again, we show in Appendix A 
that in our region E{r) is always small). The result, derived in Appendix C is: 



t(r)V3 = ty3 



2vr V rh, 



L \ /27rr 
cos 



(4.3) 



where we have assumed 

/27rro\ 



cos(^-^J=0, (4.4) 

i. e. the observer is located either at a maximum or a minimum of the density profile^^ 
and we have defined rhor(^) as the radial coordinate distance travelled by light at time 
t in an FLRW Universe (which is basically the same as in our Onion model): 

MrUt)^^-^- (4.5) 

Let us pause for a moment to understand (4.3). Firstly, we point out that the radial 
dependence of t{r) is most explicit in the expression (C.12). The reason we have 
written it in the form (4.3) is because it physically clarifies on what quantities the 
correction really depends. For example, it is nice to see how the correction to the 
FLRW result in (4.3) is controlled by the parameters of the model, e and L: e is the 
amplitude of the fiuctuations (that grows like t'^^^) and only when it becomes of 0{1) 
one can see an effect. The effect also increase as we look at higher coherence scale, 
L, although we should keep in mind that with increasing wave-lengths e decreases 
and it is really the product which is important. The largest scales that are non-linear 
today are around 60/h Mpc, so that L/r^oT ~ 0.02 and so is the product, giving 
us a correction of only a few %. We note in passing that even for scales as high as 
~ 400//i Mpc the product remains roughly the same (e ~ 0.1 and Ljr^a^ ~ 0.1) and 
therefore are equally relevant. Later in the next subsection we will see, how close to 
the observer, the "relative correction" gets enhanced from a few % to almost 20%. 



^"'^The expression could be generalized for any observer location, but it would look more cumber- 
some and it would not add anything to the understanding of the physics. 
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Figure 2: Time along the geodesic of a photon arriving at r = rp at time t = to'- the 
blue solid line is the numerical solution, the red dashed line is the analytical approximation 
of (4.3), and the black dotted line is the FLRW result. In this plot ro = 36. 5L. 

Finally, to recover the radial dependence of the right hand side in (4.3) we point 
out that the quantities e and rhor depend oxitp via (3.24) and (4.5) and the radial 
dependence oitp is given in (4.2). 

We have also solved numerically (4.1), using the solution for R{r,t) that comes 
from a numerical solution of the Einstein equations. We compare the analytical 
approximation with the numerical result in figure 2. 

We note that the agreement in the plot is very good, but in general there can 
be a mismatch when E{r) is non-negligible (for example, very far from the center). 
The same mismatch is found for z{r) and Di^r). Remarkably, however, the mis- 
match disappears when one compares two physical quantities {D^ — z), confirming 
our expectation that a linear term in E{r) should not have any physical effect. 



In the previous subsection we have found t{r) along a photon trajectory. In this 
subsection we obtain the redshift z{r) corresponding to a source located at r. The 
differential equation governing this relation is given by [3]: 



Again, ignoring E{r), we were able to find an analytic approximation. The result, 
derived in Appendix D can be written as: 



4.2 Redshift 



dz + z)B! 



(4.6) 



dr y/l + 2E 
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Figure 3: Redshift {z) along the geodesic of a photon arriving at r = ro at time t = to- the 
blue solid line is the numerical solution, the red dashed line is the analytical approximation, 
and the black dotted line is the FLRW result. In this plot rp = 36. 5L, tp = 3.3 x lO^^sec. 
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(4.7) 



Notice that the oscillating correction factor with respect to FLRW is once again 
proportional to the amplitude of the density fluctuations, e and to the scale L divided 
by the horizon size, as in (4.3). Here however the correction factor is 4 times larger, 
and also inside an exponential. 

Again we have also solved the differential equation numerically and we compare 
the result with the analytic approximation in figure 3. 

What is very interesting to notice is that for small z the correction is quite large. 
In fact expanding close to the observer we have 



vr 5r 



cos 



2Txr 
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where we have defined 



5r = r — To 



(4.8) 



(4.9) 



One can check that the pre-factor corresponds to the FLRW relation. As we had 
alluded to before, the "relative correction" near the observer is controlled by the 
ratio L/5r instead of Ljr^ar^ and this is one of the key results of the paper. Since 
near to? e ~ 1? "we see that in the first few oscillations when (5r is close to L the 
inhomogeneities provide a large 0{\) correction. This explains why several authors 
have obtained large corrections in the void models. 



21 



5 Luminosity Distance 



In this section we show how to calculate the luminosity distance in our model. As we 
said, we put the observer at a generic distance from the centre. Then we are able to 
calculate the luminosity distance as long as the source is aligned with the centre of 
the LTB model and the observer. 



5.1 Off-centre observer 

The way we proceed is to notice first that in General Relativity there is an exact 
relationship between the so-called angular diameter distance (-Da) and the luminosity 
distance (-Dl), and it is simply [44]: 

DL = DA{l + zf. (5.1) 

So we focus on the calculation of the angular diameter distance. It is well-known [3] 
that for an observer at the centre of an LTB model the angular diameter distance is 
just: 

Da = R\s. (5.2) 

where the subscript S means that the quantity has to be evaluated at the space- 
time location of the source. In order to find the distance Da for a generic location 
of the observer one has to solve the equation of motion for a bunch of geodesies 
infinitesimally close to the radial one. If the observer is not at the centre all these 
geodesies are not radial. Nonetheless we have been able to derive the exact form of 
Da in Appendix E.4, for radial sources. The final result is given by: 

Da = Rs (ro r n^9Fr ^vllVlm h ^^2 ^0 ' ^S'^) 
V Jro {l + 2E{r)){l + z{r))R{r,t{r)y J 

where the integral has to be evaluated along the radial photon trajectory. 

This expression can be easily evaluated numerically. One can get some insight 
into it by looking at the homogeneous limit (FLRW), by substituting R = a(t)r, and 
neglecting E{r). In this case the correction factor inside the brackets is: 

""Xo (l + 2E(r))(l + ^(r))/2(r,t(r))2 rs ^ ' 

This reduces to 1 in the ro — limit. Moreover, for any ro this is what we expect 
for the FLRW limit, where the distance is just 

DA = a{t){rs-ro). (5.5) 
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Figure 4: Luminosity distance {Dl) along tlie geodesic of a photon arriving at r = ro at 
time t = to- the blue solid line is the numerical solution, the red dashed line is the analytical 
approximation, and the black dotted line is the FLRW result. In this plot ro = 36.5-L, 
to = 3.3 X lO^^sec. 

5.2 Analytical Estimate 

We observe that the — z relationship can not in general be expressed as a single- 
valued function Dl{z). The physical reason is that in some region of space (collapsing 
structure with high overdensity) the photon can get a blueshift. Therefore two points 
at different distances may have the same value of z (similar to what was found in 
[45]). 

So we have a parametric expression for Dl — z, where the parameter is r. In the 
previous section we have already obtained z{r) (4.7). The expression for Diir) is 
derived in appendix (E.4): 

DLir) ^ 3/3(1 + z)5r{Mtof'^ 

We compare this analytical approximation with a numerical solution in figure 4. 

Again we comment on the size of the correction. Similar to what happens for 
z in (4.8) the correction is large close to the observer (small 5r and e of order 1). 
Remarkably the correction in DL{r) is of the same kind as in z{r) (4.8), but it is 
half its size (close to the observer). So while considering the Dl — z relationship the 
two effects will not cancel out, as one can see comparing eqs. (4.8) and (5.6) and 
comparing also the numerical results of fig. (3) and (4). The net effect is shown in the 
parametric plot of fig. (5), where we compare the analytic and the numeric solutions 
for the Dl — z relationship. 

For large redshifts, when the density contrasts are small, Dl is a single valued 



27rr 



(5.6) 
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Figure 5: Luminosity distance [Dl) vs. redshift [z): the blue solid line is the numerical 
solution, the red dashed line is the analytical approximation, the black dotted line is the 
EdS model = 1), and the green long-dashed line is the ACDM result (with Q.\ = 0.7). 
In this plot To = 36. 5L, to = 3.3 x 10^'^ sec and the observer is in an underdense point. 

function in z. For instance for large redshifts one can easily obtain Dl{z): 

1 + z — > (^-^^ ' and — > 3/3(1 + zf{Mtpfl^br (5.7) 

leading to 

Dl(^) = 3to(l + 2- v^TT^), (5.8) 

which coincides with the fiat matter dominated FLRW expression as expected, since 
at high z the density contrast is very small. 

6 Application to Cosmology and "Apparent" Ac- 
celeration 

Now, we have all the computational tools that we need for applications to cosmology. 
We stress first, that we have a cosmological model that allows to study nonlinear 
dynamics and light propagation, and this could be used to understand several phe- 
nomena (like gravitational lensing, or secondary CMB fluctuations). However our 
focus in this paper is to try and understand better the hypothesis that an inhomoge- 
neous Universe with only dust can mimick Dark Energy. 

To summarize, we have seen so far (see figs. (2,3,4)) that there is basically no 
overall effect on or on z in our Onion model: all the quantities oscillate around 
the FLRW values, with the oscillation amplitude controlled by two parameters, L/rjior 
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and e (typically a product of these two parameters). The latter is expected to become 
of 0{1) for scales in the range of tens of Mpc. However, for these scales the former 
ratio is too small, L/r^^or ~ 0.01, leading typically only to a correction of a few percent 
in the luminosity-redshift relationship. Note that a similar correction is present also 
in the evolution of the average matter density 3.8. Thus we find an overall effect 
of about a few percent coming from inhomogeneities. This is larger than what is 
computed using perturbative arguments [18, 19, 20, 21], which give a 10~^ effect, but 
still smaller than what is required to explain away dark energy. 

However, we have also seen that one cannot naively argue, based on the above, 
that inhomogeneous models will fail to give sizeable corrections (around 10%) that is 
required to explain the supernova data or any other significant effect^^. The oscilla- 
tions in the Dl — z relationship become a large effect close to the observer, because 
near the observer the relevant quantities, z and D/,, are themselves small. We saw 
that the corrections do not proportionally become smaller but instead the relative 
corrections are now approximately governed by the ratio Ljbr. Thus in the first 
few oscillation cycles when L/5r ~ OiV) we in fact notice a sizeable correction (see 
fig. (3,4)). Based on our analysis we will now discuss three main issues: 

1. How large is the effect due to a large scale fluctuation? Can it mimick Dark En- 
ergy? If it can, what is the relevant length scale L and amplitude e? 

2. Can the fact that light travels through clumpy matter give rise to any interesting 
effects? Here we really try to go beyond the Onion model, but use it to estimate the 
corrections. 

3. We discuss some observable deviations from the FLRW interpretation of the super- 
nova data coming from inhomogeneities with "typical" amplitudes (i.e. as predicted 
from CMB and/or observed in the galaxy power spectrum). 

6.1 Can a large-scale fluctuation mimick Dark Energy? 

Here we explore the possibility that a fluctuation on a (large) scale L could be re- 
sponsible for a very large effect on the Di^ — z relationship close to the observer at low 
and intermediate redshifts, such that it can mimick Dark Energy. The idea is that 
if we consider a very large scale, the local Hubble flow could be very different from 
the average one (indeed measurements at high redshifts give typically lower values 
of the Hubble constant, see the discussion in [46]). We will consider the effect of 
one single wavelength for simplicity. Then, we will vary the amplitude around (and 
above) the typical values inferred from the observed power spectrum. We will also 
vary the position of the observer with the same purpose of investigating if there is 
a setup that could mimick Dark Energy. Crucially, we note that the evidence for 
acceleration comes from the fact that we observe a mismatch between the expan- 

""^^Such a conclusion was made by [48], by studying an LTB model made of infinitely thin concentric 
shells. 
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sion at low redshift (where the Hubble parameter is measured [46], between roughly 
0.03 < z < 0.07) and the expansion at higher redshift (where supernovae are ob- 
served [47], between roughly 0.4 < ^ < 1). Indeed, we find that, depending upon 
the position of the observer (under-dense or over-dense) and the amplitude of density 
fiuctuations, it is possible to account for this mismatch. 

In the following we will discuss how to account for (i) high redshift supernova data, 
(i) local measurements of the Hubble parameter, (iii) local measurements of matter 
and baryon densities, and (iv) the position of the first acoustic peak in CMB and as 
observed recently in matter power spectrum. We end with a realistic assessment of 
this model and how one can discriminate this from the ACDM one. 



Fitting the high-z Supernovae 

For physical clarity it is useful to define an "average" model which is simply the homo- 
geneous EdS model with the same value to as the inhomogeneous one, or, equivalently, 
whose density goes as the average density in the Onion model: 

2 

(P-) = ^=P-' (6-1) 

where we denote quantities corresponding to the homogeneous model with a "bar". 
The Hubble constant today {Hq) in the average model is given by 

Ho = ^. (6.2) 

6to 

Now, the important thing to realize is that the oscillations due to inhomogeneities 
for the scales that we are interested in become negligible before we reach the realm of 
high redshift supernovae, i.e. say for z > 0.4. In this region therefore the luminosity 
distance roughly coincides with the EdS universe with a Hubble parameter Hq: 



De = ^(l + z-VT+z). (6.3) 

-no 

In order to match with observations one has to compare (6.3) with the observed 
redshift of supernovae. However, we know that, in this region, the ACDM model with 
Qx = 0.7 and /i = 0.7 is in good agreement with data. To understand analytically, 
it is therefore sufficient to compare (6.3) with Dl{z) in the "concordance" ACDM 
model given by 

= , Where Jn,{z) ^ {1 + z) f^^ , (6.4) 

So, we require that the ratio 



-Da \HqJ Jn^=Q.7{z) \h 
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be roughly equal to one. It is important to note that the ratio T> changes very little 
in the relevant range, 0.5 < 2; < 1.5: 

I?(0.5) = 0.83 >V> ©(1.5) = 0.72 . (6.6) 

Thus by choosing Hq appropriately: 

0.83 > > 0.72 , (6.7) 

we can be consistent with the high redshift supernova data. 

What is clear from (6.7) however, is that the Hubble constant of the average model 
has to be significantly smaller than what is typically measured, h ~ 0.65 — 0.72, by 
local observations, such as low redshift supernovae. In the rest of the paper, for 
analytical estimates, we will use h = 0.7 which means that to be consistent with (6.6) 
we need 

0.58 > /i > 0.50. (6.8) 

The success of the model therefore will rely on whether the local inhomogeneities can 
explain a larger local value of the Hubble parameter or not. This is what we discuss 
in the next subsection. 



Inhomogeneities & Local Hubble constant 

In order to see whether inhomogeneities can account for a larger local Hubble param- 
eter let us first compute how z is related to Dl, by combining eqs. (4.8) and (5.6). 
As we pointed out earlier, the correction in the redshift close to the observer is double 
the correction in D^, therefore we have a net effect. For small z we get 

3 / e L 27rr\ 

which is equivalent to 

z f e L 27rr\ z 

where we impose the last equality since we want to reproduce the local observations. 
Thus it is clear that, in order for our model to explain the measured Hubble parameter, 
we have to ensure, from (6.7) 

0.83 > f 1 + f ^ cos 1^ = ^ > 0.72 . (6.10) 

\ Iti or L / Hq n 

First of all we notice that one needs to consider a scale of fluctuation, L, that at least 
extends up to 2; ~ 0.07 (corresponding to a ~ l/65/iMpc~^), otherwise it is very 
difficult to ensure that the Hubble law looks linear in the range 0.03 < z < 0.08. 
This basically gives us an L > 200//i Mpc (in fact, we will see below that the optimal 
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condition is L/2 ~ 200//i Mpc). Next let us try to estimate the density contrast that 
is necessary. 

Our analytical formulas are only valid when the observer is located at a maximum 
or a minimum of the density profile, but this should be sufficient just for the purpose 
of estimating the correction^^. Let us consider these two cases separately. What we 
really want to find out is the maximal local Hubble parameter h that we can obtain 
due to the inhomogeneities, given an h. Now, we get a maximal correction (in the 
right direction) when cos(27rr/L) = —1, in which case 




(6.11) 



When the observer is at an overdense region, the smallest r for which this happens 
is 6r = 3L/4, and therefore to be consistent with (6.7) we find 

e>0.8. (6.12) 

On the other hand if the observer is at the minimal underdense region, we get a 
maximal correction at 6r = L/4 giving us a much larger correction and we only need 

e>0.27. (6.13) 

Thus, prima facie, it seems possible for the inhomogeneities to account for the ob- 
served "largeness" of the local Hubble parameter. The model prefers an underdense 
observer as one then requires smaller density contrasts and secondly for an overdense 
observer, the Hubble diagram is non-linear initially at very small redshifts z < 0.02 
which may still be conflict with observations^^. 

In fig. (6) we show one example, through numerical results, in which an almost 
linear — z relationship is recovered for z < 0.07, but a significant deviation from 
a CDM model is obtained at larger redshift, that reproduces something quite similar 
to a ACDM model. We show in fig. (7) how the customary plot m — m'^^P*^ vs. z 
looks like for these two models. In those plots m — m°™P*^ is defined as 

m - m'^P'y = 5 logio(/^L) - 5 logio(/^r"*'') , (6-14) 

where is the luminosity distance for a homogeneous FLRW model with no 

matter and no cosmological constant, but just negative curvature. The reason for 
this definition is that a homogeneous model with no negative pressure components 
will always have m — m°™P*y < 0. As one can see from the plots it is possible instead 
for an inhomogeneous model to have m — rn'^™^*^ > 0. 

Even though the model is not intended to be completely realistic (because of 
having the extra spherical symmetry and only one wavelength), we can see from fig. (7) 

^•^We could vary the position of the observer using numerical solutions, as well, but it does not 
change the situation qualitatively. 

We have not carefully studied these very local observations because local measurements of matter 
density almost certainly rules out having an overdense observer in the Onion model. 
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Figure 6: Luminosity distance (Dl) vs. redshift (z) at very small (left) and intermediate 
(right) redshifts: the blue solid line is the numerical solution of our model, the black dotted 
line is the EdS homogeneous model {Qm = 1), the green dashed line is the ACDM model 
(Qa = 0.73). The models are normalized to have the same slope at low redshift. In the 
right-hand side plot, we have superimposed the supernovae gold data set of [47]. The 
plots in the bottom show the density contrast (red line) as a function of z, for the model 
considered. 
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''^empty Q 



L = 450//iMpc ; ^/{5^ = OM 




0.25 0.5 0.75 1 1.25 1.5 1.75 



Figure 7: Magnitude residual from empty FLRW (Am) vs. redshift (z): the blue solid 
line is our numerical solution, the black dotted line is the CDM model (r^m, = 1)) and the 
green long-dashed line is the ACDM result (with = 0.73). We have superimposed the 
supernovae gold data set of [47]. 
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Figure 8: We show here the contribution to the for each data point (gold data set 
of [47]) as a function of redshift, where m is the theoretical apparent magnitude, m^xp 
is the observed magnitude and cr is the experimental error. Our model (blue boxes) is 
compared with ACDM with J^a = 0.73 (green triangles). 



Table 1: Comparison with experimental data (gold data set of [47]) 



Model 


(157 d.o.f.) 


ACDM homogeneous model (with fi^ = 0.27, I^a = 0.73) 
EdS homogeneous model (with fi^/ = 1.00, I^a = 0.00) 
Onion model {-^ (P) = 0.34 on L = 450//iMpc scale) 


178 
325 
212 



that it can provide a reasonable fit to the experimental data. We have performed also 
a analysis, and we have found the results given in tables 1,2, that show explicitly 
that the model in the example gives a result not too far from the best-fit ACDM. On 
the other hand from fig. (7) one can also see that the Onion model differs from the 
ACDM model in two regions: objects look fainter (larger Dl) at 0.1 < 2; < 0.35 and 
look brighter (smaller Dl) at z >1. With present data it is still hard to discriminate 
using this general feature: as one can see from fig. (8), individual points give similar 
contributions to the x^, when comparing our Onion model and the ACDM model. 



Table 2: Comparison with experimental data (full data set of [47]) 



Model 


X^ (186 d.o.f.) 


ACDM homogeneous model (with Qm = 0.27, ^a = 0.73) 
EdS homogeneous model (with Qm = 1.00, I^a = 0.00) 
Onion model (^(5^) = 0.34 on L = 450//iMpc scale) 


233 
403 
273 



^2 



4 

2 
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To summarize, from our study we find surprisingly that we need a density con- 
trast just slightly nonlinear {6 ^ —0.3) to account for the observed Hubble diagram. 
Nonetheless such a value represents a quite significant fluctuation about the observed 
average S that is observed on this relevant scale (|(5|~ 0.05 — 0.1). Thus it is fair to 
say that it is difficult (but perhaps possible) to explain away dark energy solely by 
including corrections coming from the first oscillation in the Onion model. Although 
the results are encouraging, probably it indicates that one has to go beyond the Onion 
model and take into account of other effects such as coming from the fact that light 
travels mostly through voids, and this is what we will discuss in the subsection 6.2. 



Matter density Si CMB acoustic peak 

Having checked that the Onion model can be consistent with supernovae data and 
the Hubble diagram, we turn our attention to other cosmological observations. 

Matter Abundance: There are different cosmological observations (for example 
the measurement of the matter-radiation equality scale looking at the matter power 
spectrum) which measure the matter density in our local universe and our model has 
to be consistent with these measurements. We already know that the ACDM model 
with h = 0.7 and Qm ~ 0.3, i.e. with 

Pmeasured = ^n^H^m% ^ 4 X lO^^^GcV" , (6.15) 

ovr 

provides a good fit to all these observations. On the other hand in our model the 
average density is given by 

3 

Po = -^Hlm^pi , (6.16) 



Q_^g ^ /Measured _ _ ^ q^^^ ^ ^g_^^^ 



SO that 

„ (~) Tj2 (-) l2 

Pmeasured ^^m-^Q \L^lb 

where in obtaining the inequalities we have used (6.7). If we identify ^measured with 
the local value of matter density, this means that it has to be approximately half of 
the average value. In particular this implies that the we (observer) must be located 
in an under dense region. 

In fact one can make an estimate of the measured local density in our model. 
Typically a measurement of this sort involves measuring the mass in a certain region 
and then measuring the volume of the same region. The crucial point to note is that 
while measuring the volume, one uses the luminosity distance as the yard-stick of 
measuring distances. Thus, for the purpose of illustration, if we consider a spherical 
region of radius 5r around the observer, then the measured density would be given 
by 

P_ed-^-^^-M,^(-j . (6.18) 
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Below we provide some numerical results for this quantity, but let us first try to 
estimate it analytically. 



Substituting (5.6) in (6.18) we find 



Pmeasured ~ r / r \ r» -i 

6n{l + zriMt,y [l-^te) 



/^measured , / r\ 27rrl3' (6.19) 



Or, approximately 

Po 

As an estimate let us consider the region upto the maximal correction (see previous 
subsection), so that cos{2nr/L) = — 1 at (5r = L/4, and we have 

(6.20) 

In order to be consistent with (6.17) we therefore find that we need 

e>0.48. (6.21) 



The example in fig. (6-7) have a local density measured inside a radius of about 100 
Mpc of l.Spm^o- We leave a more detailed investigation of checking the consistency of 
observed matter densities coming from different measurements with the expectations 
from the Onion models for future [49]. 

First acoustic peak in CMB: We have seen how our model can be consistent 
with locally observed matter density by virtue of "us" living in an underdense region. 
However, the measurements of CMB are sensitive to only the average values since 
our model reduces to the EdS model, already by a redshift of z ~ 0.3. This fact also 
ensures that our model can fit very well, at least the first peak position 

h = h{^m, hjilb/^my^-r/^m) , (6.22) 

in CMB. In (6.22) b and 7 stands for baryons and photons respectively. For instance, 
using Qm = 0.9, h = 0.58, and a standard value for the baryon-to-matter ratio, 
Clh/flm = 0.13, we find that CMBFAST [50] gives good agreement to the first peak 
position. It would be necessary to see if we can be consistent with the full CMB 
spectrum. Although this is beyond the scope of the present paper, we notice that it 
is similar in spirit to what the authors of [51] have done: by using a lower value of h 
they could fit the CMB with an EdS model. 



Baryon Oscillations: Finally our model is consistent with the recent measurement 
of the baryon acoustic peak in the galaxy distribution [52]. In fact the angular 
diameter distance at z = 0.35 looks very similar to the distance in the ACDM model 
(with Q\ = 0.7). Numerically we get for example the following values: 

L)^(0.35) = 1375 Mpc forACDM, (6.23) 
Da{0.35) = 1386 Mpc for Onion model with = 0.34 . (6.24) 



33 



It is actually very hard to compare with the measured values reported in [52]: 

Dv{0.35) = 1370 ± 64 Mpc , (6.25) 

since the "distance" Dy reported in the paper does not have a definition valid in 
general, but it is an FLRW based definition: Dy = {D\z/H{z)Y^^, where H{z) is 
the FLRW Hubble parameter. However, since our results are very close to AC DM, 
it is clear that they will be consistent also with the observations. 

Assessment 

We comment here on the viability of the model, that we have explored: 

1. As we have shown, the amplitude of the density contrast needs to be quite 
large (~ 0.3) on a scale of roughly 400/ h Mpc, with respect to the average values 
taken from the power spectrum (0.05 — 0.1). 

2. In this case, one also needs to ensure that the observer sits in a particular 
position that guarantees to have corrections of the same size, looking in different 
directions. This forces the observer to sit close to the minimum of a three dimensional 
valley (that is, close to the centre of a void). This possibility has already been 
proposed by several authors [3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13] , and we basically 
confirm the results of these previous papers, with the difference that we have a model 
that is valid from the early to the late Universe and that it is not limited to the study 
of one void. 

3. Many observations would need to be fit carefully using our parameters: for 
example one should check if a value of the matter density a little higher than the 
observed value can be consistent with all the observations. And also if the baryon 
and the matter density can be made consistent with the full CMB spectrum (as we 
said, this has been shown to be possible under some additional assumption on the 
power spectrum by [51]). 

4. It will be possible to discriminate between this Onion model and the ACDM 
concordance model with more Supernova data: in the former model objects look 
fainter (larger D^) at 0.1 ^ z ^ 0.35 and look brighter (smaller Di) at z > 1 than in 
the latter. 

Anyway the Onion model cannot be used directly as a fully realistic model, to 
interpret the data, since some crucial effects are hidden by the symmetries of the 
model. We discuss some of these effects in the next subsection. 

6.2 Light propagation through voids 

Here we try to go beyond and compensate for a couple of crucial features missing in 
our Onion model. First, in the real world the overdensities virialize; instead in LTB 
models nothing prevents the thin shells to collapse to infinite density. This means 
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we cannot follow the system beyond the time of collapse (so we can consider non- 
linear effects only at very low redshifts). Moreover, due to virialization, voids and 
overdensities behave very differently from each other in the nonlinear stage. Over- 
dense regions are prevented from collapse by angular momentum. There is nothing 
instead, that prevents a void to become more and more empty. In our model there is 
an extremely delicate balance in the Dl — z relationship between the effect given by 
voids and the effect given by structures. A photon that travels through a void gets 
significantly more redshift than the average, but this is canceled by the less redshift 
(or even blueshift!) that it gets when travelling through a collapsing region with 
high density. Now, in the real universe it is unlikely that the collapse is so strong to 
give significant blueshift. So, taking into account of this fact might well break the 
cancellation between voids and structures. Note that, if this cancellation is broken, 
this will go exactly in the direction of mimicking an accelerating Universe [49]. In 
fact, a photon travelling in the late time Universe would get more redshift than a 
photon travelling in the "linear" Universe. As a result distant objects would look less 
redshifted than close ones (as in a homogeneous ACDM model). 

The second limitation, which probably is the most relevant and interesting, is the 
fact that in the Onion model a light ray that comes from a distant source has to 
travel unavoidably through an equal number of underdense and overdense structures. 
This is not the case instead in the real Universe, where most likely a photon will 
never encounter nonlinear high density regions. In fact, only as long as the density 
fluctuations are linear, the fraction of volume in the Universe in underdensities is 
equal to the fraction in overdensities, but as they start to become nonlinear, the 
volume in voids overcomes the volume in structures. In this photon most 

likely travels through voids and therefore typically there is no compensation between 
redshift in structures and voids. Here we only present an estimate of the possible 
effect, leaving a more accurate study to [49]. 

In this context we note that in the previous literature although people have tried to 
account for changes in luminosities due to this effect, and in fact found in some cases 
significant corrections, relatively little work has been done to understand the effects 
on redshift. One of the earliest attempts dates back to Kantowski [53] who tried to 
estimate the corrections in a Swiss-Cheese model of the universe (an exact solution, 
in which spherical Schwarzschild "holes" with a mass in the centre are embedded 
in an FLRW universe), making use of the optical scalar equations [32] that regulate 
the evolution of Da ■ He estimated in this model the correction on the redshift to 
be negligible, but found important corrections coming from Da to the deceleration 
parameter (up to 50%, see also [54]). Afterwards Dyer and Roeder [29] applied the 
same optical scalar equation to a beam of light travelling in empty space (re-deriving 
the result by Zeldovich) and generalizing it to partially empty beams, confirming 
important corrections to the deceleration parameter. However, Weinberg [55] has 
shown that even if a beam of light gets a significant correction to the distance while 
traveling in empty space, another beam traveling close to a compact object gets an 
equal and opposite correction. Therefore on average the FLRW result is claimed to be 
correct, and this is interpreted as a consequence of the conservation of the number of 
photons. While this statement seems correct (the total luminosity being proportional 
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to the number of photons), it does not take into account variations on the redshift of 
the single photons (Weinberg's statement has been challenged in [45] : if the underlying 
gravitational dynamics differs from FLRW, the use of the homogeneous Dl — z is 
not justified in spite of the fact that the photon number is conserved). Weinberg's 
statement is often assumed to imply that with sufficiently many sources there in no 
need to worry about inhomogeneities. However, even accepting the statement, the 
question remains whether with a finite number of sources we may observe a result 
very different from the average. More recently, some authors [56] have studied for 
this purpose the magnification and demagnification effect using lensing techniques 
(therefore ignoring again the effect on the redshift of light due to inhomogeneous 
expansion rate) and numerical N-body simulations claiming that the dispersion in a 
realistic Universe is at most of a few percent, well below the intrinsic dispersion of 
Supernovae (and the effect can become important only at very high z). However, 
contrary to these observations in our model we have seen (both numerically and 
analytically) that the correction in z{r) is significant and must be properly taken into 
account while deriving the luminosity vs. redshift relation. 

In fact the correction in z{r) is approximately double that of in Diir) and im- 
portantly in the opposite direction. So in our attempt to try to estimate the effect 
of "photons travelling mostly through voids" in Dl{z), we are only going to try to 
estimate the correction to the redshift in a given cycle, had it not passed through an 
overdense region, and then divide it by half. As one can see from a z{r) plot (figure 9), 
as the photon approaches a highly overdense structure, it first suffers from a blue shift 
(let's call the coordinate range through which this happens Ar) and then takes some 
distance (approximately equal to Ar) to recover from this blue shift. Thus as an 
estimate we see that the photon receives a "zero-shift" in a coordinate range 2Ar. 
Thus it is natural to conjecture that if the photon had not encountered the structure 
it would not have lost the 2Ar coordinate distance and therefore the increment in the 
redshift over a cycle should be corrected as 

{Az),orr ~ (A2;)o„ion (^1 + = (A2;)plrw (^1 + . (6.26) 

The second equality follows from the fact that in the Onion model over a period, i. e. 
length L, the change in the redshift is the same as the homogeneous EdS model. 

We note that this effect on the redshift (applied to CMB photons) has been esti- 
mated already by Zeldovich and Sazhin [57] in the context of a Swiss-cheese model. 

Let us see whether we can estimate the crucial ratio Ar/L. It is clear that the 
blue-shift begins at a maximum of z{r) and it ends at the minimum of the curve. 
Starting from the expression of z{r) in (4.7) one finds 



ar 



27rr\ L ( lixr 
2 sm —— H cos 



, ~ y J > 1- (6-27) 

The second term in the right hand side is much suppressed because of the ratio L/r^on 
and therefore we approximately have 

f2Txr\ 1 

sm — =-. (6.28) 
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Figure 9: Redshift(z) along the geodesic of a photon arriving at r = ro at time t = to in 
a strongly nonlinear regime: the blue dotted line is the numerical solution, the red dashed 
line is the analytical approximation (4.7), and the black solid line is the FLRW result. In 
this plot ro = 36L, to = 3.3 x lO^^sec. 



This has periodic solutions signalling maxima and minima with an interval of 




(6.29) 



As we mentioned before, the correction is really expected to be half of the above due 
to the cancelling effect in Dilr) and therefore we find 

Ar 1 1 _i / 1 

Correction ~ — = sm — 

L 2 TT V2e 

For e = 1, we get a maximal correction of 1/3 ~ 0(1), while as we go back in 
time, the density contrast decreases and so does the correction. In fact, we find no 
blue-shift at all once e < 0.5, which occurs approximately at the redshift 

l + z = — ^ zr^l, (6.31) 

precisely what is required to explain the high-redshift supernovae! 

In general if we consider a small scale, that becomes nonlinear earlier, this already 
can give some effect at higher z. 



(6.30) 



6.3 Minimal uncertainties in the — z plot 

In this section we discuss what are the minimal implications of our results on the 
distance-redshift measurements. From what we have seen so far, it is perhaps possible 
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that the effects are large enough to mimiclc Dark: Energy, but it can also be that dark 
energy really is the explanation for the Di — z relation. Even in this case, it is of 
great importance to study what are the deviations from FLRW. 

We find two main effects: first, a generic observer will see a correction (at all 
redshifts) to — z of roughly 0.15 apparent magnitudes with respect to the naive 
FLRW result, second any observer will see more scatter in the low redshift supernovae 
than in the high redshift ones. 

To show this, we just take realistic values of the density contrast 6 on some spe- 
cific scale. For simplicity we just use the same setup previously described, but we 
decrease the amplitude to the measured average value. To be more precise one would 
need to include fluctuations on all scales, but we leave this for a subsequent paper 
[49]. So, we consider here an amplitude ^/{S^ ~ 0.075 (where (...) is an average 
over many domains) on a scale L 400//iMpc, that is roughly consistent with the 
measured matter power spectrum. We choose one particular position of the observer 
(underdense) and we compute numerically the deviations from the FLRW case. We 
find uncertainties of order 0.15 apparent magnitudes'^, which is the same size as the 
uncertainty on the intrinsic luminosity of a single Supernova ([47]). Crucially, how- 
ever, while the uncertainty in the intrinsic luminosity can be decreased by increasing 
the number of measured supernovae (the error decreases as 1/V^), our result does 
not decrease since it is a real physical effect that depends on the mass distribution 
around us and it is therefore the most relevant systematic effect (for a discussion of 
the uncertainties see for example [58]). 

We show the "Hubble" diagram of the example analyzed in fig. (10) and its mag- 
nitude compared to an empty universe in fig. (11). We show the "error" Am that 
is done using a CDM model with respect to this example in fig. (12) (we expect a 
similar result even in the presence of a dark energy component). 

We note also that such an uncertainty is of the same order of the detected 
anisotropy of the Hubble constant at low z [59]. It is usually assumed that the 
anisotropy is due to the velocity of our galaxy with respect to the FLRW frame 
(although with controversial results [59]), but our results imply that the large scale 
structure plays a role that may be responsible for at least part of the anisotropy. 
Similarly one could study the effect on the CMB due to the large scale structure (see 
for example [57], [60]). 

Finally the other fact that we can account for is the experimental observation that 
low redshift supernovae show more scatter in apparent magnitudes than high redshift 
one [31] (see also [15]). As we have seen, in fact, the corrections in z and in Dl are 
larger close to the observer. We show as an example a Di — z plot with a smaller 
wavelength {O{50)/hMpc) in fig. 13. 

-'^^A correction of this order was noted by [15], but it was not recognized that the corrections can 
affect all measurements, not just the low-redshift ones. 
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Figure 10: Luminosity distance (Dl) vs. redshift {z) at very small (left) and intermediate 
(right) redshifts: the blue solid line is the numerical solution of our model, the black dotted 
line is the EdS homogeneous model (f^m = !)> the green dashed line is the ACDM model 
(r^A = 0.73). The models are normalized to have the same slope at low redshift. The plots 
in the bottom show the density contrast (red dashed line) as a function of z, for the model 
considered. 
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L = 400//iMpc ; ./{5^ = 0.075 
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Figure 11: Magnitude residual from empty FLRW (Am) vs. redshift {z): the blue solid 
line is the numerical solution, the black dotted line is the CDM model (0^ = 1)) the green 
long-dashed line is the ACDM result (with i}\ = 0.73). We have superimposed the red thin 
dashed line, that shows the density contrast seen by the photon along its trajectory as a 
function of z. 
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Figure 12: The black solid line shows the correction in apparent magnitudes due to a 
realistic inhomogeneous density fluctuation {VP ~ 0.075, for L = AOO/hMpc) with respect 
to a homogeneous FLRW model. We have superimposed the red thin dashed line, that 
shows the density contrast seen by the photon along its trajectory as a function of z. 
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Figure 13: The blue solid line shows the oscillations in the Dl — z relationship with a 54//i 
Mpc nonlinear density fluctuation. The black dotted line is the EdS model with the same 
high-z behaviour. The green-dashed line, as usual, is the best-fit ACDM model {VL\ = 0.7). 

7 Conclusions 

In this paper we have proposed an exact model, based on LTB metric with only dust, 
for structure formation applicable to a Universe homogeneous and isotropic on very 
large scale. The structures here look like concentric shells, and we could study their 
evolution from the early Universe to the present epoch. 

In this background we have studied exactly the propagation of light along a radial 
trajectory, and we have computed the distance-redshift relation as measured by an 
observer sitting in a generic position. 

For this purpose we have derived (for the first time) an expression for the luminos- 
ity distance in an LTB metric for an off-centre observer, we have solved the system 
numerically, and we have found accurate analytical approximations to the dynamical 
equations and to the light propagation equations. We have made several conclusions 
on the basis of our computations. 

First we notice that the corrections become larger close to the observer (far away 
the corrections due to underdensities tend to cancel the corrections due to overdensi- 
ties) . 

Then, we have investigated whether there is a setup that can mimick an accel- 
erating ACDM cosmology. We have shown that this is possible under some special 
conditions: the observer has to be located around a minimum of the density contrast, 
and has to live in a large-scale underdensity with quite high density contrast: a typ- 
ical example is 5 ~ 0.3 on a scale of roughly 400//iMpc, while the typical value of 6 
at that scale is around S ~ 0.05 — 0.1. 
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Without such a large density fluctuation the situation is still very interesting: one 
can estimate a typical effect, which affects all measurements, on the Dl — z plot of 
about Am ~ 0.15 apparent magnitudes. By considering fluctuations of smaller length 
scales, ~ 100//i Mpc, we can also explain the observation that low redshift supernovae 
should show more scatter than high redshift ones (as seen experimentally [31]). 

Finally, and most crucially, we notice that making the model more realistic could 
uncover relevant effects: the fact that light in a radial trajectory in the Onion model 
has to go through an equal number of structures and voids leads, in fact, to peculiar 
cancellations in the — z plot. Without these cancellations (which do not happen 
in the real world) we have estimated that the effect goes in the direction of apparent 
acceleration and that it should be large when matter fluctuations become nonlinear. 
We leave a more detailed analysis of this effect to a subsequent paper [49] . 
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A Validity of small u and small E approximation 

The analytical analysis in our paper relied on the u variable being small. This is 
always true if one focuses on "large" enough values of r. Let us now make this 
statement more precise. 

In our model E[r) is given by (3.21) 
Now, since E[r) is a periodic function, we have from (3.3) 



9v^V,^„)V3,/M. (A.2) 



TT / V 7rr 



Thus the small u approximation breaks down approximately at r = r„(t) given by 



2 , X 1 

3 / n fF:\ 3 



ru{t) = ^ p- (Mt)3 ^ ^ = ^(Mt)3 . (A.3) 

7r\7r/ L \ n n 
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From (A. 3) it is clear that is a monotonically increasing function of time and since 
we are interested in describing the evolution only till today, we have 



2 



TT / vr 



In other words, as long as r > r„ „iax, the small u approximation is valid. We can 
actually estimate r^ max more precisely by substituting Ai from (3.28): 



2 

L ira \tcj IT \ tj ^ ' 

Thus as long as our observer is located at a position rp > QAL, we can trust our 
small-M approximation. In the numerical simulations we actually place the observer 
at a much larger distance, tq ~ SOL. 

Let us now check when E{r) remains small. Clearly this function grows with r 
and therefore we have 



M\ro + r^o.)L 



TT 



(A.6) 



where r^or is the coordinate interval that a CMB photon covers in its journey towards 
us. Now, we can estimate this number using the flat FLRW model (it is expected to 
change only slightly in the LTB model that we are considering). We find 

M..„.(^. (A.7) 
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Thus we have 



2 2 

3 



Or, 



where we have defined 



^ = —. (A.9) 

Since we are considering inhomogeneity scales roughly corresponding to 400 
Mpc//i, while the horizon is around 3Gpc//i, the ratio^^ Ljr^ax ~ 0.1 is not very 
small. In fact neglecting £'(r) introduces some mismatch between the analytical and 
numerical results in Dj^ir) and in z{r\ However, as we have discussed already in 
section 4.1, the mismatch disappears when comparing physical observable quantities 
(such Dl vs. z). For smaller scales, instead, there is no mismatch at all. 
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For large distances the "proper distance" R oc r. 
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B Near- Centre Region 



We have focused on E{r) of the form (3.21) because we want to describe the large scale 
periodic inhomogeneous structures. However, such an ansatz is actually inconsistent 
with the analytic requirements of E{r) near r = 0, and in particular in [17] it was 
claimed that this requirement does not allow for "accelerating" cosmologies. Since 
we put our observer at large distance away from the centre, it is perhaps obvious 
that this result is irrelevant for us, but here we provide one particular consistent E{r) 
which for large r reduces to that of the Onion model. 

It was pointed out in [17] that one requires the curvature to have the following 
form near r = 

E(r) ~ r" with n > 2 . (B.l) 

A simple function which has this property but asymptotes to (3.21) for large r is 
given by 

^f--) = (;^rr^ j (^ J (^) . (B.2) 

By construction 

J A.M^rL gi^2 for r > 6 

E{r) — I ^„ (^4^) sill' if) for r < 6 " ^^"^^ 

We note in passing that near the centre u tends to increase (u — »• oo also cor- 
responds to the curvature dominated regime). Thus we feel it is easy to attribute 
special properties to the central point and one of the reasons why we wanted to avoid 
placing the observer at the centre. 

It is possible to avoid u becoming large by, for instance choosing 

''>^(^V. (B.4) 



but this also tends to wash-out the inhomogeneities. 

C The photon trajectory: t(r) 

We want to solve the evolution equation for the photon trajectories given by (4.1). 
As usual we ignore E{r) in (4.1). Then, substituting R' we have for large r 

^ = -3/?(Mt)i [l + a(Mt)lA(r)\ . (C.l) 
dr L J 

It is convenient to work with the dimensionless variable 

x={Mty/^. (C.2) 



44 



The differential equation then reads 



dx{r) 
dr 



Or, 



dx = -dr (3M 



1 + aAix'^ sin 



1 + aAxX^ sin 



27rr 

~17 



2nr 



(C.3) 



(C.4) 



One cannot solve analytically this equation but we will adopt an adiabatic ap- 
proximation scheme where we treat the function x{r) (or equivalently t{r)) to be 
slowly varying with respect to the oscillatory functions. Thus while differentiating 
or integrating, whenever we encounter a product of x{r) (or one it's powers) with an 
oscillatory piece, we always neglect the variation of x{r) and treat it as a constant. 
The physical reason why this becomes a good approximation is because while the 
variation in x{r) is determined by the Hubble scale (inversely proportional), a vari- 
ation in the oscillatory functions is inversely proportional to the coherence scale L. 
Thus as long as L/rhor is small we are safe. 



Applying this method to the differential equation (C.4) we find 



x{r) 

where we have assumed 



{Mtoy/^ - /3M{r 



ro) 



aAiL 



27i(r - ro) 



-x^ cos 



27rr 



cos 







(C.5) 



(C.6) 



i.e. the observer is located either at a maximum or minimum of the density profile. 
First note, that we recover the FLRW limit by letting Ai 0. In this case we find 



x{r) = xo — [3M{r — ro) = xpir) 



(C.7) 



which is indeed the relation for an FLRW universe. Let us now check explicitly that 
(C.5) is indeed a good approximate solution to the differential equation (C.4). 



Differentiating (C.5) with respect to r we obtain 

dx{r) 



dr 



Ai — (ML)xcos 
vr 







(T)] 





aAix^ sin 



27rr 



Comparing this with (C.4) we find that (C.5) should be a good approximate solution 
provided the extra second term (relative to (C.3) is small: 



a3 ~ 

w = Ai — {ML)x < 1 , 
vr 

Since, this is again a monotonic function in time we have 

(3{ML){Mto)^/^ 



(C.^ 



Sr. 



{Mt, 
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where we have further used (3.28). In order to estimate this quantity note that the 
horizon at a given x can be estimated using just the FLRW resuh (C.7). Using (3.29) 
we therefore find 



5r 



1 / L 



TT V'"lior(a;o)y \tc J TT V^hor(xo) 

Thus, (C.5) is indeed a good approximate solution. 
Solving (C.5) we find 



to 



< 1 



(C.9) 



x(r) 



2/5a^cos(^) 



so that 



(C.IO) 



(C.ll) 



We have checked that (C.ll) is in excellent agreement with the numerical solution. 
In fact keeping upto linear order in Ai in (C.IO) 



x(r) ^ xf{t) 



1 + oiA\XF\r] cos 



(67r 



(C.12) 



we already get an excellent approximation. 

It is instructive to write the above expression as 



x\r) ^ xpir) 



e I L \ ( 2Txr 
1 + — I 1 cos . ^ 



2vr V '^h, 



which makes the dependence of the correction on the parameters of the model appar- 
ent. 



D Calculating the Redshift 

In this appendix we try and obtain the redshift z{r) corresponding to a source located 
at r. The differential equation governing this relation is given by 

dz_ ^ {i + z)ii 

dr V1 + 2E ■ 

Again, ignoring E{r), we have 
dz 



l + z 



R'dr = 2PM [x~^{r) + 2ax{r)A{r)] dr 
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2/5M 



MLaAi 2nr , . . 

■ sin — — h 2ax{r)A[r) 



{Qtt 



dr , 



where we have used (C.12). The first term can be integrated in a straightforward 
manner to yield the FLRW result. The second term can also be integrated easily. For 
the third integration one can use the same trick of first integrating with respect to 
r assuming that x is a constant and then put back the r-dependence of x. After all 
these integrations we find 

ln(l + z) -\nC = -2\n XF{r) 



{MLfaAi (2Tir 
sin 



2 ~ 27rr 

—afiA\MLx{r) cos — i— , 

Lj 



(D.2) 



37r(67r)i/'^ " \ L J n 

where C is an integration constant. Now, the last two terms are only important close 
to ts ^ to- During this time the second term is always a lot smaller than the third 
because ML <^ x» and thus we have 



ln(l + z) -\nC = -21na;F(r) 



2a/5 , , , 27rr 

-AiMLxyr) cos — — 



TT 



We can fix the integration constant by demanding that at r = ro, z = 0. This gives 
us 

lnC = 21nxo, (D.3) 

so that we have 



ln(l 



-2 In 



XF{r) 
xo 



2a(3 . ~ ^ , , 27rr 
-AiM Lx{r) cos 



TT 



(D.4) 



Or, 



1 + zir) 



X. 



o 



exp 



-AiM LxFyrj cos 



X, 



o 



TT 

2ap 



TT 



AiMLxf(t) cos 



L 

2'Kr 



(D.5) 



Again, it is instructive to write the above expression as 



l + z 



X 



o 



Xp[r) 



exp 



2e f L 



COS 



2'Kr 



X, 



o 



Xp{r) [ TT \rii, 



L 



cos 



'27rr ' 
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E Dl for an off-centre observer 

In this appendix we derive an expression for the luminosity (or angular) distance 
as seen by an observer that looks in a radial direction and that sits in a generic 
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off-centre position, at coordinate tq- In order to do tliat, we will need to consider 
trajectories that end up at the observer space-time position: however we need not 
only radial trajectories, but also those that have an infinitesimal deviation angle. 
Then, given a solid angle 6Q at the observer for a bunch of light rays that arrive 
from the source, we will be able to compute the cross sectional area A along the past 
light cone. This gives us the "angular diameter distance" (or the "area distance") 
Da = X' As we said, this is related immediately to the luminosity distance of a 
source by Da = (1 + z)'^Dl. Finally we will double-check our derivation by using the 
optical scalar equation derived by Sachs in 1961 [32, 53] in appendix F. 



E.l Quasi-radial trajectories 

First of all, we give the geodesic equations that we will use to study the trajectories 
of the photons towards the observer. The relevant geodesic equations are given by 
[24]: 

dH R'R' fdr^^ 



,o + +RRC-' = 0, (E.l) 

dv^ l + 2E\dvJ ^ ' 

d'^e ^^R! drde ^^Rdtde ^^^q ^^^q ('^^^ _ g (E 2) 

dv'^ R dv dv Rdv dv \dv J ' 

^ + 2— — — + 2- — — + 2cote— — = (E3) 
dv"^ R dv dv R dv dv dv dv 

where v is an affine parameter, and where we have dropped the equation for r{v), 
since we will not need it. In addition, the condition {ds^ = 0) gives us: 



dtV R'^ fdr^^ 



dv 1 + 2E \dv 



R'C = , (E.4) 



where 



'^^-(ij -'"'nil ■ 



Now, first we choose the observer to sit in the 9 = tt/2 plane. And then we observe 
that for symmetry reasons a trajectory lying in this plane will never escape it: we 
can in fact consistently set 9 = n/2 in the geodesic equations as all the terms in (E.2) 
vanish. The axial symmetry of the system also ensures that it is sufficient to consider 
only such trajectories for the purpose of computing the luminosity distance. So, (E.3) 
then simplifies to 

d'^cp ^ ^R' dr d(j) ^ ^R dt d(f) ^ (E 6) 

dv'^ R dv dv R dv dv 
Again, without loss of generality we can assume the observer to lie along the ray 
= 0. The radial trajectory will then be characterized by having (f){v) = 0. 



48 



However, in order to compute the luminosity distance, we need to find another 
"near-by" trajectory which, ahhough reaches (r = ro, = 0) at t = to-, starts, say at 
t = ti, with = (pi. The parameter (pi, in fact, parameterizes a bunch of geodesies 
which reach the same destination but diverge in the past. Now, we will only need 
to find trajectories which are infinitesimally close to the radial geodesic, i.e. (pi <^ 1. 
This provides a great simplification, as we notice that ~ (pf and therefore r(A, (pi) 
and t{v,(pi) do not get any correction ~ 0{(pi). To see this, observe that in general 
for any coordinate 



(E.7) 



For (p{v,(pi) this means 

(p[V,(pi) = (pi — 



SO that 



dv J ' \dvd(pi 



(Ei 



<f>i=Oy 



Thus the equations (E.l) and (E.8) do not get any correction of the order ~ (pi and 
therefore 

r{v, 0,) = r{v, 0) + O(0^) and t{v, 0^) = t{v, 0) + 0(0^) . (E.9) 



E.2 OfF-centre observational coordinates 

Before calculating the luminosity distance we need to do an additional step: we 
define the system of coordinates used by the observer (the so-called "observational 
coordinates" [38]). They are defined by the following conditions: 

I. The time-like coordinate is the proper time of the observer 

II. The first spatial coordinate is a monotonic parameter along the past light cone: 

it can be an affine parameter (v) or the coordinate r — or anything else. 

III. The remaining two coordinates are "angular coordinates" {0,(p), in the sense 
that the metric along the light cone (keeping fixed the other two coordinates) 
looks like: 

ds^ = v\d9^ + d4>'^sme^) . (E.IO) 

For the time coordinate, we can directly use t, as it is already the proper time. 
Then we attack the problem of finding the correct angular variables by first defining a 
rigid translation of the spatial coordinates, centering them in the observer, and then 
rescaling the angles in such a way that (E.IO) is satisfied. 

Let us first concentrate on a coordinate transformation from the centre of the LTB 
patch, C, to the observer point located somewhere off the centre, O. The original 
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LTB coordinates are called = {t, r, ^, 0}. We assume the observer to be comoving 
with the dust. In reality our Galaxy may have a peculiar velocity relative to the local 
LTB flow, but for the sake of simplicity we ignore it. That is, we assume the observer 
to be at a fixed coordinate distance tq from the centre. Now, let us define a new 
coordinate system with O as the origin, denoted by = {t,r,6,(j)}. The following 
relations hold between these two set of coordinates: 

t = t, 

r sin 6 sin = f sin 6 sin , 

r sin 6 cos = ro + r sin 6 cos , 

rcos6' = fcos6'. (E-H) 

Without loss of generality we assume the three points C, O, and S to build our plane 
of ^ = 6* = |. The zero reference line of and in this plane is then defined by the 
hue C-0. 

Our goal is now to determine the metric in the new coordinates, using 

The non-vanishing partial derivatives that we will need to calculate the area distance 
follows: 

dr 

dQ f 
^ 1 0=0= 1,^=0 = - , 



90, 



= f,</'=0 ^ 



(E.13) 



The angular part of the metric at the observations point looks like: 



^"2 



ds\r=const = R'^^idO' + dcj)' siu 9') . (E.14) 



At this point, in order to satisfy (E.IO) we redefine coordinates by: 

t = i, 

f = f , 

ir0 = 

Ke = e + {K-i)^. (E.15) 

In this way the angular part of the metric looks like (around 6 = 7r/2): 

-2 

ds\,=,onst = K^R^y^{dd^ + rf0') . (E.16) 
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So we require 



K = Um^^oQ) = ^, (E.17) 

dv 



where we have used the fact that [24, 38]: 

| = -(1 + .). (E.18) 

In this way the barred coordinates are observational coordinates. 

Finally, we will make use of the transformation between angles in the two coordi- 
nates systems: 



arccos 



ro + r cos{K(()) 



y^VQ + r"^ + 2ror cos{K(j)) 
where the last equality holds for infinitesimal angle. 



^ Kef), (E.19) 



r + ro 



E.3 Calculating the luminosity distance 

By definition, the area distance is the ratio between the solid angle in the observational 
coordinate and the area along the past light-cone [38]: 

2 _ dA desd(ps . .-P 

''- = ^ = Jo^v^^ (E.20) 

where as usual 5* stands for source and O stands for observer, and where we have 
already chosen Q = j1. This can be rewritten as: 

0\='^R\s- (E.21) 

Wod<po 

We have therefore to compute the ratios 

dds d<^s 



d9: 



o a<Po 



(E.22) 



for infinitesimal angle around = (which corresponds to aligned Centre, Source 
and Observer). By axial symmetry it is sufficient to consider one ratio: the two ratios 
are the same. We will therefore consider the ratio ~- . 

d(po 

In order to find what this ratio is, we consider the geodesic equation (E.6). 
Remarkably it can be simply written as: 

So, the following quantity is a constant along the geodesic: 

R'^ = C. (E.24) 

dv 
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So we may integrate the equation as: 

dv 



s 
o 



s = C I (E.25) 



where we have put = at the observer. We basically already have the result, since 
we have the final angle and so the final area A. Explicitly, we may change variable of 
integration along the geodesic, instead of the affine parameter v. Using (E.18), we 
can express the integral as: 

The last thing to be worked out is the constant C. Making use of (E.19) we have 
that: 

d(f) d(f) df ^ d(j) d(j) d(j) dr dt ^ d(f) dcj) ^ 
dv df dv d(j) dv df dt dv dcj) dv 

In the last equality we have used the fact that at the first order it is sufficient to 
consider a radial trajectory, where df = dr. Explicitly evaluating the derivatives at 
f = gives: 



Q d(j) dr dt 
Q ' df dt dv 



o ' 



(E.28) 



which means that C has the value: 

C = Ro^o . (E.29) 

This allows us to calculate the ratio of (E.22): 

d6s n R'irMr)) 

d<l)o Jro {l + z{r))y/l + 2E{r)R{r,t{r)y 

So the final result luminosity distance Dl = (1 + z^D, is given by: 



D, = il + zfRs I Ro r ,^ , , 1 • (E-31) 

J ro 



[\ + z{r))^\^-2E{r)R{r, t[r)f 



With a very good approximation we can ignore Eir) with respect to 1: 

I)^ = (1 + zfRsRo { r 7TT^7Wr^nV2dA = (1 + zfRsRoI • (E.32) 
\Jro [1 + z{r))R{r,t{r)y J 



E.4 Analytical Approximation 

In this section we will simplify the expression for luminosity distance (E.32) for the 
Onion model. In order to do this first we will evaluate approximately the integral /. 
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Using (3.22), (3.23) and (D.5) we find 



R! 



^xpAiML cos (^) + aAixj, sin (^) 



(E.33) 



It is easy to see that for the realistic L, the third term always dominates over the 
second and therefore can be ignored. Moreover while integrating the third term 
one can once again use the adiabatic approximation {{xf/tY is changing slowly as 
compared to the sine), so that we straight-forwardly get 



5r 



aAix%Lro ( ^irr 
1 : z cos 



2'Kr6r 



L 



(E.34) 



Substituting I in (E.31) we get 



DL = a + zfRsRo 



Sr 



SPrroXQ 



aAiXpLro ( 27rr 
1 ; — cos 



27rr5r 



Using the expression for R (3.5) in the small u approximation we find 

r A. rr2 

Dl ^ 3P{1 + zY6rx] 
^ 3/3(1 + z)Srxl 



aAix%Lro ( 27ir 
1 z cos 



4 



2'7ir6r 
aAiXpLro 



27ir5r 



cos 



L 

2'Kr 



(E.35) 



where in the last approximation we have used (D.5). This gives us the luminosity 
distance as a function of r. Coupled with (D.5) we can obtain a parametric plot of 
Diir) vs. z{r). 

Finally, noting that the factor in the denominator is the same quantity that ap- 
pears in the density profile we can rewrite Dl as 



Dl = 3/3(1 + z)6rx 
^ 3/3(1 + z)Srxl 



eLro ( 2'Kr 
cos 



2'Krbr 



L 

2Txr 

~17 



(E.36) 



where the last approximation is only valid for the source close to the observer, but 
typically this is when the corrections are important. 



F Derivation using the Optical Scalar Equation 

In this subsection we show that our result (E.31) is consistent with the so-called 
optical scalar equations derived first by Sachs [32]. This equation has been used for 
example in the literature to study the distance Da when light travels through empty 
regions [28, 53, 29]. 
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The Sachs equation for the evolution of a cross sectional area A for a cone of light 
originating from one point is 

y^ = 7^ where 7^ = . . (F.37) 

Here / = a/A, Rfj,^, is the Ricci tensor and is the derivative of the position of the 
photon along its path: 

^^,^dx^^^^ dtjv) dr{v) ^ \ 
dv \ dv ' dv ' ' / 

We want to see if our result obeys this equation, in order to have an independent 
double-check. 

We may express the derivatives in k^^ using again (E.18) and (4.1), so that: 



fc^= -(l + ^),^^^^^^(l + ^),0,0 . (F.38) 



It is clear from the definitions that the quantity / is proportional to the area 
distance Da- Therefore Da has to obey the same equation (F.37). 

Now, it is straightforward to check that our solution (E.31) obeys the following 
equation: 

^fDA^^fR 

Da dv^ Rdv^' ^ • ^ 

So we have to compare this with the term TZ of (F.37); by direct computation the 
relevant components of the Ricci tensor read: 

Ru = 2^ + ^, 
" R R'' 

1 R'(2E'(r) - 2RR' - R'R) 

Rrr = TTTT— — F.40 

l + 2E{r) R ^ ' 

So, contracting with k^ it is easy to see that in fact TZ = "^J^ . 
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